Two-pulse interference and superluminality 
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Abstract 

We examine how the interference of a coherent hght-pulse with its shghtly time- delayed copy 
may generate a pulse nearly identical to the original one and ahead of it. The simplicity of this 
2-pulse system enabled us to obtain exact analytic expressions of the pulse distortion, valid for a 
wide class of pulse shapes. Explicit results are given for the pulses usually considered (gaussian, 
hyperbolic secant) but also for more realistic pulses of strictly limited duration. We finally show 
that the efficiency of the 2-pulse system is comparable to that of the other superluminal systems, 
at least for the pulse advancements actually demonstrated in the optical experiments. 
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Although unambiguous demonstrations of the phenomenon have been achieved in the 
1980's Q], the now called superluminal propagation^of light pulses (fast light) still attracts 



considerable interest. For reviews see, e.g., [3|, U, O]. In an ideal experiment, the pulse 
having covered some distance L in a linear medium with anomalous dispersion leaves it 
with almost no change in shape, in significant advance of a pulse which would have covered 
the same distance L at the velocity c of light in vacuum. As a matter of fact, the field 
transmitted at some time t is not the consequence of the signal entering the medium at the 
same time t but that of all the signals anterior to t by more than L/c and this puzzling 
phenomenon is not at odds with the relativistic causality. It is now well recognised that the 
superluminal propagation originates in a predominantly destructive interference between the 
different frequency-components of the pulse, all propagating with velocities smaller than c 
(see, e.g.,j6|). A similar mechanism occurs in systems involving a non-dispersive medium 
when several modes or paths of light-transmission, with different propagation times, are 
possible. The interference of the different transmitted pulses may then generate an output 
pulse ahead of the fastest one among them and even in absolute advance of the input 
pulse [3, 18|]. This explains in particular the superluminal transmission of optical pulses 
through a multilayer structure j^, Q though, due the high multiplicity of possible paths, 
the phenomenon is usually analysed in terms of evanescent waves (optical tunnelling). A 
simpler and more direct demonstration of multiple-pulse interference is provided by the 
recent experiment involving a fibre taper coupled with a microsphere HI]. The output 
pulse then results from the interference of the pulses transmitted by the taper after having 
achieved n roundtrips (0 < n < cxo ) inside the sphere (whispering-gallery mode). In order 
to examine how a time-delay can originate a time- advancement, we examine in this paper 
the still simpler situation in which only two pulses interfere (2-pulse systems). This occurs 
in particular in birefringent media where the light can propagate at two different velocities. 
Superluminality has then been inferred from phase measurements by Solly et al. [12l] and 
directly demonstrated by Brunner et al. jisl ]. 

The causality principle fixes a strong link between the transmission-delay of a linear 
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system and its gain profile [5|, 114, 1151, 116| . The largest superluminal effects are generally ex- 
pected when the frequency uoq of the optical carrier coincides with a well-marked minimum 
of gain. In the present case this is actually achieved when the 2-pulse interference is predom- 
inantly destructive, as indicated before. The relation between the slowly-varying envelopes 
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of the incident and output pulses ej„(t) and eoutif) reads then Coutit) = Aein(t) — Bemit — r) 
with A > B > {] and r > 0. The coefficient A (B) is the amphtude-gain for the fast 
(slow) pulse and r is the delay of the slowest pulse with respect to the fastest one {t is 
thus a local time). Finally the 2-pulse system can be characterised by its impulse re- 
sponse h{t) = A6{t) — B5{t — t) or, in the frequency domain, by its transfer function 
H{Q) = A — B exp{—iQT) , where Q is the deviation of the current optical frequency u from 
cuq. Our definitions and sign conventions are those used in signal theory [17]. The phase 
shift = aTg[H{Q)] and the amplitude-gain |-ff(f2)| are respectively odd and even pe- 

riodic functions of Q , the latter oscillating between m = A — B (its value for i7 = 0) and 
M = A + B (constructive interference). FigH] shows the frequency dependence of |if(0)| 
and of (p{^l) in a case taken as reference in the following. 
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FIG. 1: Amplitude-gain |i?(J7)| (full line) and phase shift (p{0,) (dashed line) of the 2-pulse system 
as a function of the detuning Q = uj — ujq , for M = 7.4 and m = 0.89. The frequency unit is 1/r. 
The spectrum of a gaussian pulse of duration Tp = 13r is given for reference (dotted line). 



Some light on the superluminality mechanism is obtained by doing a Taylor series expan 
sion of Coutit). At the ffist order in r, we get 

dCin it) 



Coutit) ^{A- B)eUt) + i-B)i-T)^^ ^ m 



ein{t) + ai- 



dt 



(1) 



with ai = ~ l] f • This leads at the ffist order in ai to Coutif) ~ mcinit + ai) and 
H{yi) K, mexp(if2ai). Despite their apparent roughness, the previous approximations are 
exactly those currently used to define the group velocity, ai is nothing but the group time- 
advancement dip/dVL\Q^=Q and is actually positive since the negative signs associated with 
the time-delay r of the slow pulse and with the destructive interference compensate one 
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another (see the intermediate form in Eq{T]). At this order of approximation, the output 
pulse has the same shape that the input pulse, is m times larger and is ahead of the fastest 
transmitted pulse by ai. The approximation will be good if mexp{iQai) is itself a good 
approximation of the exact transfer function on the entirety of the pulse spectrum. This 
requires (at least) the time-delay r to be small compared to the duration Tp of the pulse 
(hereafter the full width at half maximum of its intensity profile). The examination of the 
expression of ai then shows that group advancements comparable to the pulse duration can 
be attained only if the maximum gain M is large compared to the gain m at the frepuency 
Uo of the optical carrier. This property is common to all the superluminal systems [18']. In 
purely passive (not amplifying) optical systems M < 1 and the condition M/m » 1 is 
only met when m << 1. A large post-amplification (often not explicitly mentioned) is then 
required to normalise the amplitude of the output pulse. Conversely some post-attenuation 



may be required in active systems to obtain m ^ 1 [19|. In both cases the device achieving 
the normalisation of the output pulse should be considered as a part of the superluminal 
system and the gains considered in the following (as well as in FigH]) are those of the overall 
system. Note however that the double condition M >> m ^ 1 can be met in all-optical 
systems, for example in an amplifying medium with a dip in its gain-profile, as proposed in 



Iq and demonstrated in |20|. 



Superluminal effects are always accompanied by some pulse distortion. The distortion 
is negligible only when the spectrum (duration) of the input pulse is very narrow (long) 
but the fractional advancement is then also negligible. Although this point is often not 
addressed, the real challenge is to attain significant fractional advancements with moder- 
ate pulse distortion. The latter is conveniently characterised by the root-mean-square (rms) 
deviation D of Coutit) from the ideal advanced form ej„(t + a) , where a is the aimed advance- 
ment. By normalising the energy of the input pulse {J^^ \ein(t)f dt = 1), D is reduced to 

1 /2 

D = Icoutit) — ein{t + 0.)!^ dt^ . Minimising the distortion is in general a formidable 

problem without analytical solution ISj. The 2-pulse system considered here is probably 
the only system for which exact analytical results can be obtained. For sake of simplicity, 
we assume that emit) is real (no chirping). Coutif) is then also real and D depends on the 



pulse shape only via p(r 
autocorrelation of ej„(t) 



j_p(a) and p{a + r), where p{t) = emit + 9)ein{0)d9 is the 



13]. Recall that p{t) is an even function of t, has its maximum 



(equal to 1) for t = 0, and is the inverse Fourier transform of the energy spectrum of ej„(t). 
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For given delay r and advancement a, simple calculations show that D attains a minimum 



, p2(a) + p2(a + r) - 2p(r)p(a)p(a + r) 
(a, r ) = ^ / 1 l-p2(r) 

when M = and m = ^^^^^j^^j^y^- The distortion only cancels in the trivial case 

a = 0. We have then M = m = 1 , that is A = 1 and 5 = (no interference). For a > 0, 
the minimum distortion is attained for m < 1 and the group advancement ai generally 
differs both from a and from the advancement of the pulse maximum. It keeps however 
a physical meaning. A simple application of the so-called moment theorem in signal analysis 



171] indeed shows that ai is the advancement of the centre-of-gravity of the pulse envelope, as 

I ! 

large as the distortion may be [16] . If ein{t) is continuous, the derivatives p'^^^t) = d"'p/dt'^ 
of p(t) are continuous up to n = 2 and the distortion has an absolute minimum Damin when 
r ^ 0. It reads 

Da^^n {a) = + (3) 

where p*^^^(0) < 0. At this limit, m = p(a), ai = p(-^|p(2"(o) and M = t^Ist^^- The latter 
expression shows that, for a given advancement, the lowest distortion is attained when the 
maximum gain M of the system tends to infinity. When M is large but finite, we easily get 
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WW) 



-p{a) 



(4) 



Mp(2)(0) 

that is Dminici, M) / Damin{,o) = 1 + 0(1/M). The asymptotic limit Damin can thus be ap- 
proached for moderate values of M. The above results are valid under the unique assumption 
that the pulse envelope emit) is continuous. Since discontinuities in dem/dt may originate 



strong transients in the output pulse (see, e.g.. Fig. 3 in 21|), we will now consider the usual 
case of bell-shaped envelopes, the derivative of which is continuous. The autocorrelation p(t) 
is then also bell-shaped, its derivatives p*-"-* (t) are continuous up to n = 4 and it is possible 
to obtain approximate expressions of Damin, ^ and ai/a valid at the 2nd order in a. They 
readD„™„(a) ^ ^p(4)(0) - [p(2)(0)]', m ^ l+p^^\Q)^ and ^ ^ 1 + f _p{2)(o) 
At this level of approximation, the distortion is mainly due to the frequency-dependence of 



the gain (Fig HI), w. 
of the pulse itself 



rich originates a broadening of the pulse spectrum and thus a narrowing 
Due to the pulse-area-conservation theorem |22], this narrowing is 



accompanied by a pulse magnification and this explains why the rms distortion is minimised 
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when the gain m at Q = (that is at the centre of the pulse spectrum) is shghtly smaller 
than unity. 

The previous results apply to a wide class of pulse shapes. We first consider the 
gaussian pulses, the most popular ones in the literature on superluminality. We write 
their envelope under the form ein{t) oc exp (-tV2), such that Tp = 2V\n2 ^ 1.67 
and p{t) = exp(— We then get Daminia) = a/1 + (1 + exp (— 0^/2) and and 
Oi = a. Note however that the common value of ai and a differs from the advance- 
ment ttm of the pulse maximum. The latter is easily determined from the expression 
of eout{t) and reads a^, = ^a/1 + 4af — ij /2ai. At the 2nd order in a, we finally get 
Daminia) ~ 0^/2^2 ~ 0.98 {a/Tpf with m ^ 1 — a? /A and am/a ^ 1 — a^. 

As a second example we consider the classical hyperbolic-secant pulses (sech-pulses) 
ein{t) oc 1/cosht, such that Tp = 2 argcosh(v^) ^ 1.76 and p{t) = t/sinht. We then 



get Damin{ci) ~ (sir^) ^ (coth a — ^) +1 and ai = 3 (cotha — 1/a) 7^ a , that is, 

in the 2nd order approximation, Damin ~ 2a^/3-\/5 ~ 0.93 {a/Tp)^ , with m ~ 1 — a^/6 and 
am/a ^ 1 — a^/15. 

Contrary to the analytical pulses above considered, the pulses actually used in the ex- 
periments have obviously a strictly finite duration and their envelope has thus at least two 
points of non-analyticity. There is then an unavoidable distortion due to the fact that, in 
our local time picture, the output signal cannot anticipate the beginning and the end of the 
input signal. Moreover, when M is very large, these points of non-analyticity may generate 
strong transients in the output pulse 2l|], reinforcing the distortion. All that occurs even if 
the envelope of the input pulse is quasi-gaussian. In the sequel we will take as reference the 
"cos-pulses" of envelope ej„(t) oc (1 + cost), for \t\ < n and ej„(t) = elsewhere (total width 
27r). They are such that Tp = 2arccos(V2 - 1) ^ 2.29, p{t) = (l - S) + %r 

for |t| < 2tt and p{t) = elsewhere. We now get Daminict) = a/1 — p^(a) — 3/2(a), 
ai/a = 3f{a)/ap{a) ^ 1 and = arctan(ai), with /(a) = (l - ^) ^ + • At 

the 2nd order in a, these expressions read Damin ~ a'^/3y/2 ^ 1.23 (a/rp)^, ai/a ~ 1 and 
am/a ^ 1 — a^/3, with m ^ 1 — a^/6 . We observe that, for moderate fractional advancements 
a/Tp, the distortion of the cos-pulses in the 2-pulse system is not much larger than that of 
the standard gaussian pulses. Note that, in both cases, the advancement is smaller than 
the group advancement. 

As indicated before (see EqJl]), the asymptotic limit Dam,in can be approached for reason- 
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able values of the maximum amplitude-gain M. For further comparisons, we consider the 
case a/Tp = 0.32 and M = 7.4. The corresponding gain and phase profiles are those of FiglH 
The distortion remains moderate (14%), exceeding its absolute minimum Damin only by a 
factor 1.25. We compare FigI2]the envelope Coutit) to its ideal form Cinit + a). As expected. 




FIG. 2: Envelopes ein{t) of the input cos-pulse and eoutit) of the output pulse in the 2-pulse 
system for a/Tp = 0.32 and M = 7.4 (full line). The ideally advanced envelope ein{t + a) is given 
for comparison (dotted line). The time unit is the duration Tp of the input pulse. 



one can see large deviations at the beginning and the end of the input pulse. Occurring 
in the pulse pedestal, these features are obviously much less apparent on the correspond- 
ing intensity-profile (FigO]), which is the profile usually observed in the optical experiments. 
Apart from a slight overshot at the end, the profile of the output pulse is comparable to that 
obtained with an input gaussian pulse of same duration. In both cases, the main distortion 
consists in a narrowing of the pulse, resulting in different advancements on its fall and on 
its rise. 

In the abundant literature on superluminality, one finds few direct demonstrations of 



significant pulse-advancements (say a > Tp/10) with moderate distortion 
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26 



27 



28 



23 



2j, 



29| . Indeed, problems of noise, of instability and of hypersensitivity to small 
defects in the input pulse severely limit the maximum gain M which can be actually used 
in the experiments (keep in mind that we include in the system gain that of the device 
eventuahy used to normalise the output pulse, so that m ^ 1 in the moderate distortion 
limit). In our experiment at a wavelength A ~ 3 mm j^, we succeeded in using M 160. In 
;he optical domain, quantum effects (amplified spontaneous emission, detection noise, etc) 
3o| further reduce the gains actually usable and, to our knowledge, all the experiments have 
been achieved with M < 10 (amplitude gain). Due to these restrictions on the maximum 
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FIG. 3: Normalised intensity-profiles corresponding to the pulse envelopes ej„(t) and eout{t) of 
Figl2] (full line). The profiles obtained with a gaussian pulse of same duration Tp are given for 
reference (dashed line). 

gain, the efficiency of a superluminal system may be defined as its ability, M being given, 
to attain a fractional advancement a/r^ as large as possible with a moderate distortion. The 
acceptable distortion obviously depends on the advancement, very small when the fractional 
advancement is itself very small 3l| and the larger the larger a/r^ is. In the following we 
retain as acceptable the distortion observed in our reference case (Figl2]and FigJSj), namely 
D 0.45a/rp . FigH] shows the fractional advancements a/r^ attained in this way for various 
values of M ranging from 1 (no advancement) to 20. The advancement actually measured 
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FIG. 4: Fractional advancement a/rp attainable with an acceptable distortion (see text) as a 
function of the maximum gain M (including the gain of the normalisation device) . The input pulse 
is a cos-pulse. The full (dotted) line and the full (open) circles relate to the 2-pulse (single-dip) 
system. 



in the experiments is the advancement a™ of the pulse maximum (a^ < o) and the most 
accessible parameter is the contrast Tmax/To between the maximum intensity-transmission 
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of the optical device and its value at ujq. Since Tmax oc and Tq oc m^, we obviously have 
Tmax/To = M"^ /rn? ^ M^. The dependence of am/Tp on Tmax/To is shown Fig|5l 




FIG. 5: Fractional advancement am/Tp of the pulse maximum as a function of the ratio Tmax/To 
of the maximum intensity-transmission of the optical system over its value at the centre of the 
pulse spectrum. The input pulse is a cos-pulse. The full (dotted) line and the full (open) circles 
relate to the 2-pulse (single-dip) system. 



For comparison consider now the sup erluminal systems havirig demonstrated the most 



striking fractional advancements 
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19|, 



23 



24 



25 



26 



271, l28|, l29|. Most involve an 



homogeneous medium, the transmission curve of which has a dip at the frequency cjn. 



dip may be natural 



0,Q 



but is often induced by an additional light field [25 



26 



271, 



28 



'he 



29| 



(which unfortunately may also induce extra noise). Using our general conventions, H{Q) 



Mexp 



i+in/-f 



is generally a good approximation of the transfer function of these single- 



dip systems. In this expression Z is the field optical-thickness associated with the dip 
and 7 is the width of the corresponding line, assimilated to a Lorentzian The group 
advancement, the amplitude gain at = and the intensity-contrast Tmax/To respectively 
read ai = Z/^y, m = Mexp(— Z) and Tmax/To = exp(2Z). The distortion can easily be 
expressed in the frequency domain by means of the Parseval's theorem [l^. We get [16] 
D = ^J^^\H{n)-exp (ina) | ^ | (^]) | ^ c/fi/27r where ^^„(^^) = e,„(t) exp(-ifit)c/t is 
the Fourier transform of ej„(t). M and a being fixed, the parameters Z and 7 are adjusted in 
order to minimise the distortion. This is efficiently achieved by means of a genetic algorithm 



32l |. Rapid convergence is obtained by exploring values of Z and 7 close to those giving 
m = 1 and Oi = a, that is Z = In M and 7 = In M/ a. For each value of M, the fractional 
advancement a/Tp is calculated by a dichotomy method: (i) Choose < a/Tp < 1 (ii) 
Minimise D by the genetic algorithm (iii) If Dmin is too large (small), decrease (increase) 



a/Tp. Go to (i). We stop the calculation when Dmin equals the acceptable distortion within 
0.1% and we determine the corresponding advancement am of the pulse maximum. FigH] 
(FigE]) shows a/Tp (am/Tp) as a function of M (Tmax/To) [33|. Comparing the advancements 
ttm attained with the 2-pulse system and with the single-dip system, we see that the former is 
only 22% smaller than the latter for the largest ratio Tmax/To actually involved in an optical 
experiment (Tmax/To ~ 100). In most experiments Tmax/To < 10 and the two advancements 
are nearly equal. Conversely the single-dip system would be much more efficient if ratios 
Tmax/To >> 100 could be used. Quite similar results are obtained by considering the other 
superluminal systems having demonstrated significant fractional advancements, in particular 
the periodic fibre Bragg gratings j23], the taper coupled with a microsphere 11 1 and the 



doublet of gain lines 



16|, 



19| . A better efficiency is expected from a doublet of absorption lines 



34| but, to our knowledge, this system has not been actually used. 



To summarise, we have studied the superluminal system involving the predominantly 
destructive interference of a pulse with a time-delayed copy. Due to the simphcity of this 
system, we have succeeded in obtaining exact analytical expressions of the distortion. Finally 
we have shown that, despite its simplicity, this system permit to observe superluminal effects 
comparable to all those which have been actually demonstrated in optics. 

Laboratoire PhLAM is Unite Mixte de Recherche de I'Universite de Lille I et du CNRS 
(UMR 8523). CERLA is Federation de Recherche du CNRS (FR 2416). 
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Insofar as the prevailing distortion originates in the 2nd order variations of the gain, Dmin oc 
/ In M for a gaussian pulse 16|, so that a oc {Dmin/o) In ^ ■ The nearly linear curve of FigH] 



(dashed line) shows that this is approximately valid for the cos-pulse in the range 1 < M < 20 



Using again a genetic algorithm, we find that for our reference gain M = 7.4 the fractional 
advancement of the pulse maximum would be am/Tp = 0.44 instead of 0.32 with the single-dip 
system and 0.29 with our 2-pulse system. 



12 



